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ABSTRACT 
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The  optical  method  of  caustics  as  applied  to  the  measurement  of  the  J  integral  in  plasti¬ 
cally  deforming  materials  was  studied  numerically  and  experimentally.  The  numerical  cal¬ 
culations  modeled  a  cracked  body  loaded  under  Mode  1,  plane  stress,  small  scale  yielding 
conditions.  The  results  were  used  to  generate  simulated  optical  patterns  for  caustics 
obtained  from  both  inside  and  outside  the  crack  tip  plastic  zone.  These  caustics  were  com¬ 
pared  with  corresponding  caustics  observed  experimentally.  In  addition,  the  region  of 
dominance  of  the  plane  stress,  HRR  field  was  investigated,  and  the  effect  of  plasticity  on 
the  accuracy  of  caustics  from  the  elastic  region  outside  the  plastic  zone  was  determined. 
The  extent  of  the  region  of  three  dimensionality  near  the  crack  tip  and  its  effect  on  the 
method  of  caustics  was  also  investigated. 
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1.  INTRODUCTION 

In  recent  years  the  optical  method  of  caustics  has  been  developed  into  a  very  success¬ 
ful  experimental  tool  for  studying  linear  elastic  fracture  mechanics  problems.  In  such 
problems,  the  method  of  caustics  can  be  used  for  the  direct  measurement  of  the  stress 
intensity  factor,  K,  [1,2].  Particularly  important  applications  of  caustics  have  been  made  in 
dynamic  fracture  mechanics  where  the  stress  intensity  factor  cannot  generally  be  deter¬ 
mined  analytically.  The  past  success  of  caustics  has  led  to  efforts  to  extend  the  method  to 
applications  in  elastic  plastic  fracture. 

Some  preliminary  work  [3,4],  based  on  the  assumption  of  the  validity  of  the  plane 
stress,  HRR,  asymptotic  crack  tip  field  [5,6],  demonstrated  that  the  value  of  the  J  integral 
can  be  directly  measured  with  caustics.  Unfortunately  the  region  of  dominance  of  the 
plane  stress  HRR  field  has  not  been  accurately  established.  Thus  the  conditions  under 
which  the  analytical  results  reported  in  [3,4]  are  valid  are  uncertain.  Experimental  results 

given  in  [4,7,8]  indicate  that  the  method  is  promising  and  is  worthy  of  further  investiga¬ 
tion. 
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In  this  paper  the  method  of  caustics  is  investigated  experimentally  and  numerically. 
One  of  the  goals  of  this  work  is  to  provide  an  analysis  of  caustics  which  is  based  on  a  full 
field,  small  scale  yielding,  finite  element  solution.  In  small  scale  yielding,  well  outside  the 
plastic  zone,  the  elastic,  singular  stress  field  dominates.  Inside  the  plastic  zone,  very  near 
the  crack  tip,  the  HRR  field  dominates.  In  the  transition  region  between  these  two  fields  no 
analytical  solution  is  known.  Thus  no  analysis  of  the  method  of  caustics  exists  for  this 
region.  This  places  certain  limitations  on  the  applicability  of  caustics.  To  quantify  the  limi¬ 
tations,  the  extent  of  dominance  of  the  plane  stress  HRR  field  is  studied  numerically.  In 
addition,  an  analysis  of  caustics  based  on  the  numerical  results  is  introduced.  This 
analysis  is  not  limited  by  the  assumption  of  a  particular  asymptotic  field. 

For  caustics  obtained  from  the  elastic  region  surrounding  the  crack  tip  plastic  zone,  an 
analysis  based  on  the  linear  elastic  Kt  field  may  be  used.  However  crack  tip  plasticity  often 
affects  caustics,  causing  errors  in  the  measurement  of  K,.  The  resulting  errors  are  investi¬ 
gated  here  experimentally. 

The  analysis  of  the  method  of  caustics  applies  to  cracked,  planar  bodies  that  are  thin 
enough  so  that  the  state  of  deformation  is  one  of  plane  stress.  However,  there  is  a  region 
near  the  crack  tip  where  the  deformation  fields  are  three  dimensional.  Experiments  were 
performed  to  investigate  the  extent  of  the  region  of  three  dimensional  fields  for  plastically 
deforming  materials,  and  to  determine  the  effect  of  three  dimensionality  on  the  accuracy 
of  the  method  of  caustics  as  based  on  plane  stress  fields. 


2.  THE  METHOD  OF  CAUSTICS 

Consider  a  set  of  parallel  light  rays  normally  incident  on  a  planar,  reflective  specimen. 
The  specimen  surface  was  initially  optically  flat  but  it  is  deformed  due  to  the  tensile  crack 
tip  loading.  The  deformed  shape  of  the  specimen  surface  is  such  that  the  virtual  extension 
of  the  reflected  light  rays  forms  an  envelope  in  space  as  illustrated  in  Tig.  1.  This  surface, 
called  the  "caustic  surface"  is  the  locus  of  points  of  maximum  luminosity.  Its  intersection 
with  a  plane  located  a  distance  z0  behind  the  specimen  is  called  the  "caustic  curve".  This 
curve  bounds  a  dark  region,  called  the  "shadow  spot".  By  placing  a  camera  in  front  of  the 
specimen  to  collect  the  reflected  light  rays  and  by  focusing  the  camera  behind  the 
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specimen  the  caustic  can  be  photographed. 


Suppose  that  light  reflected  from  a  point  (x.ltx2)  on  the  specimen  intersects  the  plane  at 
z0  at  a  point  (X1,X2),  where  (x!,x2)  is  a  coordinate  system  on  the  specimen  centered  on  the 
crack  tip  and  (XlfX2)  is  a  system  translated  by  a  distance  z0  behind  the  specimen.  Then 
{XltX2)  are  given  by  the  mapping  [1] 


Xn  xa  2z0 


df(x  i,X2) 


dxn 


Q  =  1,2, 


(1) 


where  f(xr,x2)  =  -i/3(xi,x2),  the  out  of  plane  surface  displacement  of  the  specimen.  As  dis¬ 
cussed  in  [9]  a  caustic  curve  will  exist  if  and  only  if  the  Jacobian  determinant  of  the  map¬ 
ping  vanishes, 


J(xhx2,z0)  =  det 


dx3 


=  0. 


(2) 


The  locus  of  points  on  the  specimen  satisfying  J  =  0  is  called  the  "initial  curve".  All 
points  on  the  initial  curve  map  onto  the  caustic  curve.  In  addition  all  points  inside  and 
outside  the  initial  curve  map  outside  the  caustic.  Since  the  light  that  forms  the  caustic 
curve  is  reflected  from  the  initial  curve,  whatever  information  is  conveyed  by  the  caustic 
comes  from  this  curve  only.  Equation  (2),  defining  the  initial  curve  depends  on  z0.  Thus 
by  varying  z0  the  initial  curve  position  may  be  varied.  If  z0  is  large  then  the  initial  curve 
will  be  located  far  away  from  the  crack  tip.  If  z0  is  small  then  the  initial  curve  will  be  close 
to  the  crack  tip. 


When  the  initial  curve  is  well  outside  the  crack  tip  plastic  zone  the  analysis  of  caustics 
based  on  the  linear  elastic  K,  field  may  be  used.  In  such  cases  the  Mode-1,  plane  stress,  out 
of  plane  displacement  field  is  given  by 


u3 


-uhKj 

EA27it 


cos 


2  ' 


(3) 


where  E  is  the  elastic  modulus,  u  is  Poisson’s  ratio,  h  is  the  specimen  thickness,  and  r  and  9 
are  polar  coordinates.  By  substituting  Eq  (3)  into  Eqs  (1)  and  (2)  it  is  found  that  the  caustic 
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is  an  epicycloid  and  that  K,  is  related  to  the  caustic  diameter  D  (width  of  the  caustic  in  the 
X2  direction)  by  [1,2] 


ED5/2 
10 .7z0vh  ' 


(4) 


The  initial  curve  is  circular  and  its  radius  r0  is 


r0  =  0.3 16D 


3hi>KiZ0 

2V2  ;rE 


2/5 


(5) 


When  the  initial  curve  is  within  the  region  of  dominance  of  the  plane  stress  HRR  field  an 
alternate  analysis  of  caustics  must  be  used.  For  a  material  with  a  stress  strain  curve  that 
may  be  described  by  the  Ramberg  Osgood  model 


where  a0  and  c0  are  the  yield  stress  and  yield  strain  in  tension,  the  asymptotic,  out  of 
plane  displacement  is  given  by  [5,6] 


oa0h 

JE 

n 

K+l 

2  E 

oallnr 

Err{0,n)  +  Egs(0,n) 


(6) 


where  /„  is  a  numerical  factor  that  depends  on  the  hardening  exponent  n  and  Err  and  E0l) 
are  dimensionless  functions  of  9  and  n.  Substitution  of  Eq  (6)  into  Eqs  (1)  and  (2)  yields  [3] 


J  =  Sn 


ovj 

E 


E 

oa0  za  h 


fl  +  1 

n 


3n+2 

D  "  , 


(7) 


where  S„  is  a  numerical  factor  dependent  on  n.  For  a  non-hardening  material  [4] 


(8) 


-  6  - 
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The  caustic  for  a  hardening  exponent  of  n=9  is  shown  in  Fig.  2.  Figure  2a  shows  the 
predicted  caustic  and  Fig.  2b  shows  the  caustic  observed  in  an  experiment.  The  plastic 
zone  can  be  seen  surrounding  the  caustic  of  Fig.  2b  demonstrating  that  the  caustic  came 
from  light  reflected  from  within  the  plastic  zone. 

Unlike  the  elastic  case  the  initial  curve,  sketched  in  Fig.  3,  is  no  longer  circular;  it 
depends  on  the  hardening  level  of  the  material.  The  point  on  the  initial  curve  that  maps  to 
the  maximum  value  of  X2  on  the  caustic  curve  is  located  at  an  angle  9max  from  the  axis. 
Let  the  distance  to  this  point  from  the  crack  tip  be  called  r0,  the  initial  curve  size.  It  is 
important  to  know  the  location  of  this  point  since  it  is  mapped  to  the  maximum  diameter 
of  the  caustic,  which  is  used  in  the  interpretation  of  the  results.  It  was  shown  in  [8]  that 
$max  varies  from  72°  to  56°  as  n  varies  from  1  to  oo.  For  a  power  law  hardening  material  [8] 


r„ 


f 

0.385 D  n= 9 
0.40D  n= 50 


(9) 


Note  that  Eqs  (3)-(9)  are  based  on  the  assumption  of  the  validity  of  particular  asymp¬ 
totic  fields.  In  the  next  section  an  attempt  is  made  to  eliminate  this  restriction  by  con¬ 
structing  caustics  based  on  a  full  field  numerical  solution. 


3.  NUMERICAL  ANALYSIS 

3.1.  Formulation. 

The  numerical  calculations  modeled  a  semi-infinite  crack  under  Mode-I,  plane  stress, 
small  scale  yielding  conditions.  The  displacements  from  the  singular,  elastic  crack  tip  field 
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ujr.0)  =  K, 


ua(0) 


were  specified  on  a  circle  of  radius  of  approximately  3400  times  the  the  smallest  element 
length,  as  shown  in  Fig.  4.  The  maximum  extent  of  the  plastic  zone  was  contained  within 
1/30  of  this  radius.  The  cutout  in  Fig.  4a  is  a  fine  mesh  region  around  the  crack  tip  which 
is  shown  in  detail  in  Fig.  4b.  An  incremental  J2  plasticity  theory  was  used.  The  material 
obeyed  the  von  Mises  yield  criterion  and  followed  a  piecewise  power  hardening  law  in 
uniaxial  tension  of  the  form 

a  <a0 

n 

<?  >  ao 


with  hardening  exponent  n=9.  This  value  was  chosen  to  match  the  4340  steel  used  in  the 
experiments.  All  plasticity  was  confined  to  the  active  region  (Fig.  4a)  consisting  of  1704 
four-noded  isoparametric  quadrilaterals.  The  quadrilaterals  were  formed  from  4  constant 
strain  triangles  with  static  condensation  of  the  internal  node.  Static  condensation  was 
used  in  the  region  around  the  active  mesh  that  always  remained  elastic.  The  stress  compu¬ 
tations  were  performed  using  a  tangential  stiffness-radial  return  method  with  subincre¬ 
mentation  [10].  The  nonlinear  finite  element  equilibrium  equations  were  solved  incremen¬ 
tally  using  the  iterative  Newton-Raphson  method  [1 1], 

After  the  element  nearest  to  the  crack  tip  has  yielded,  the  load,  applied  through  the 
stress  intensity  factor  Kh  was  increased  in  steps  of  5-10%  of  the  value  at  yield  until  the 
extent  of  the  plastic  zone  ahead  of  the  crack  lip  was  50  limes  tire  smallest  element  size. 


3.2.  Results. 
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The  calculated  and  experimentally  observed  plane  stress,  crack  tip  plastic  zones  are 
shown  in  Fig.  5.  A  point  was  plotted  for  every  element  that  had  yielded,  thus  defining  the 
calculated  plastic  zone  in  Fig.  5a.  Figure  5b  shows  visual  evidence  of  plastic  deformation 
observed  on  the  surface  of  a  thin,  compact  tension  specimen  of  4340  steel  with  a  harden¬ 
ing  exponent  of  n=8.7.  It  is  seen  that  the  calculated  and  observed  plastic  zones  agree  in 
shape.  Both  calculations  and  experiment  show  that  the  extent  of  the  plastic  zone  ahead  of 
the  crack  tip  is  rp=0.2S(K,/ao)z. 

The  stress  distribution  oZ2/t0  ahead  of  the  crack  is  shown  in  Fig.  6.  The  non- 
dimensionalizatlon  of  by  (K,/a0)2  is  used  because  the  plastic  zone  scales  with  this 
parameter.  The  computed  stress  agrees  with  the  stresses  of  the  plane  stress  F1RR  field  to 
within  1%  when  x^O.08 {K,/a0)z,  i.e.  x^O-Sr,,.  In  Fig.  7  the  plastic  strain  distributions 
622 Ao  and  <h/(0  are  compared  to  the  strains  predicted  by  the  F1RR  solution.  Again  it  is 
seen  that  the  numerical  results  agree  with  the  HRR  feld  up  to  xt~0 3rp.  Of  particular 
interest  is  the  result  for  f33/(0  since  this  is  related  to  the  out  of  plane  displacement,  which 
is  important  to  the  method  of  caustics. 

The  stresses  ahead  of  the  crack  are  compared  to  the  linear  elastic  field  in  Fig.  8.  At 
the  boundary  of  the  plastic  zone  (xi=0.25(K,/<to)2)  <t22  begins  to  decrease  rapidly,  making  a 
transition  to  the  linear  elastic  solution  as  x,  increases.  It  is  seen  that  the  linear  elastic 
field  agrees  with  the  actual  stress  field  when  xt  >  1.5  rp. 

3.3.  Numerically  Generated  Caustics. 

To  provide  a  means  to  analyze  caustics  that  is  not  dependent  on  the  assumption  of 
dominance  of  the  HRR  or  K,  fields,  caustics  were  generated  using  the  results  of  the  FEM 
analysis.  The  out  of  plane  surface  displacements,  shown  here  in  Fig.  9,  were  smoothed 
using  a  least  squares  scheme  [12]  and  caustics  were  generated  by  mapping  light  rays  point 
by  point  using  Eq  (1)  for  different  values  of  z0.  The  caustics  are  shown  in  Fig.  10  for 
values  of  rQ/rp  from  0.19  to  1.3.  It  is  seen  that  for  r0/rp= 0.19  the  numerically  simulated 
caustic  agrees  in  shape  with  the  caustic  predicted  using  the  the  HRR  field,  Fig.  2.  When 
rQ/rp=  1.3  the  numerically  simulated  caustic.  Fig.  lOf,  agrees  with  the  caustic  predicted 
using  the  elastic,  field.  In  the  region  between  r0/rp=.  19  and  1.3  there  is  a  transition 
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from  the  "HRR  like"  caustic  to  the  "elastic  caustic",  it  is  seen  that  for  r0/rp  as  small  as  0.30 
(Fig.  10b)  the  caustic  shape  deviates  from  the  HRR  shape. 

The  numerical  caustics  were  produced  for  a  fixed  level  of  applied  K,,  or  J  since  J=K?/E. 
By  measuring  the  caustics  a  relationship  between  the  caustic  diameter  and  the  known 
value  of  J  was  determined  for  all  values  of  r0/rp.  The  relationship  is  given  in  Fig.  1 1  where 
the  non-dimensionalized  caustic  diameter  is  plotted  as  a  function  of  the  non- 
dimensionalized  J  integral.  The  bars  on  the  numerical  results  indicate  the  uncertainty  in  D 
due  to  the  discretization  of  the  finite  elements.  The  solid  line  gives  the  relationship  for 
caustics  from  the  K,  dominant  region,  Eq  (4)  with  0.3.  The  dashed  line  gives  the  rela¬ 
tionship  for  caustics  from  the  HRR  dominant  region,  Eq  (7),  with  ck=  1.0.  The  value  of  a=1.0 
was  chosen  to  match  the  material  idealization  of  the  numerical  model.  Note,  however,  that 
Eq  (7),  is  a  very  weak  function  of  o-. 

For  large  values  of  z0  (large  initial  curves)  the  numerical  results  approach  the  elastic 
relation,  Eq  (5).  For  small  values  of  z0  (initial  curve  close  to  the  crack  tip)  the  numerical 
results  agree  with  the  relation  obtained  on  the  basis  of  the  HRR  field,  Eq  (7). 


4.  DESCRIPTION  OF  EXPERIMENTS 

To  experimentally  determine  the  range  of  validity  of  the  caustics  measurements,  the  J 
integral  (Jeans)  or  the  stress  intensity  factor  ( Kcaus )  was  measured  with  caustics  using 
different  initial  curves.  These  quantities  were  simultaneously  measured  from  the  boundary 
conditions  and  are  denoted  by  JBC  and  KBC.  The  values  were  then  compared  for  different 
conditions. 

The  caustics  were  photographed  using  the  apparatus  illustrated  in  Fig.  12.  The  system 
consisting  of  the  lens  and  screen  forms  a  telephoto  camera  that  is  focused  on  a  virtual 
image  plane  at  a  distance  z0  behind  the  specimen.  Light  reflected  from  the  specimen 
forms  a  caustic  on  the  screen  which  is  then  photographed  with  a  35  mm  camera.  Graph 
paper  on  the  screen  allows  the  magnification  from  the  screen  to  the  film  to  be  determined 
and  the  thin  lens  equation  is  used  to  calculate  the  magnification  from  the  virtual  image 
plane  to  the  screen.  By  varying  the  distance  from  the  lens  to  the  screen  and  by  using 
different  focal  length  lenses,  the  system  can  be  focused  at  different  distances  z0  behind 


-  10  - 


the  specimen.  Varying  z0  varies  the  initial  curve,  allowing  data  to  be  obtained  from 
different  distances  from  the  crack  tip. 

Note  that  to  have  large  z0  one  must  make  the  distance  from  the  lens  to  the  screen  small. 
In  order  to  maintain  a  good  magnification,  a  long  focal  length  lens  should  be  used  when  z0 
is  large. 

The  initial  curve  was  calculated  from  the  caustic  diameter  using  Eq  (5)  for  r0>rp  and  Eq 
(9)  for  rQ<rp.  In  transition  regions  where  neither  the  HRR  or  the  elastic  field  dominates,  or 
in  three-dimensional  regions,  these  equations  are  approximate.  However,  r0  should  always 
be  bounded  by  the  two  values,  0.3 16D  and  0.40D,  (Eqs  (5)  and  (9))  and  thus  any  calculated 
value  of  r0  using  these  equations  is  approximately  correct  even  in  the  transition  regions. 

For  these  tests,  loading  beyond  relatively  small  scale  yielding  conditions  was  not 
applied.  Thus  only  the  applied  load  needed  to  be  measured  to  determine  J.  Values  of  KBC 
(and  hence  JBC)  were  calculated  using  formulas  given  in  [13]  of  the  form 

K[  =  f  {P, a', geometry)  (10) 


where  P  is  the  load  applied  to  the  specimen.  A  correction  to  the  crack  length  to  account 
for  the  plastic  zone  in  calculating  K,  was  used.  The  correction,  discussed  in  [14]  is  given 
by 


a'  =  a  +  <j>ry 


where  a  is  the  crack  length, 


1  n- 1  ,Ki 


Ky  2k  n+l 


and 


[1  ~{P/Plm)2Y 


- 11  - 


Equations  for  the  limit  load,  PLM,  are  given  in  [14]. 

The  specimen  dimensions  and  geometries  are  given  in  Fig.  13.  The  test  specimens  were 
lapped  optically  flat  and  polished  to  a  mirror  finish.  The  pre-crack  was  cut  with  a  wire 
electric  discharge  machine,  producing  a  crack  tip  diameter  of  0.3  mm.  Annealed  4340  steel 
and  cold-rolled  1018  steel  were  used.  The  heat  treatment  of  the  4340  steel  was:  675°C  for 
1.5  hours,  oil  quench,  then  anneal  at  520°C  for  1  hour  and  air  cool.  The  material  proper¬ 
ties  were:  yield  stress  a0  =  830  MPa,  hardening  exponent  n=8.7,  and  a  =  3.3  for  a  best  fit 
of  the  uniaxial  tension  properties  to  the  Ramberg-Osgood  material  model.  The  1018  steel 
was  used  as  received.  It  had  a0  =  560  MPa  with  no  hardening. 

The  tests  proceeded  by  loading  the  specimen  to  some  level  and  then  maintaining  that 
level  while  caustics  were  photographed  for  different  values  of  z0.  The  load  was  then 
increased  a  step  and  the  process  repeated. 


5.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 

5.1.  Sequence  of  Caustics. 

A  sequence  of  photographs  of  caustics  is  shown  in  Fig.  14.  These  photographs  are  from 
specimen  28,  loaded  such  that  rp=3.4mm.  Only  the  distance  z0  was  varied,  thus  varying 
the  initial  curve  size  r0.  The  parameter  r0/rp  in  the  figure  is  the  ratio  of  the  initial  curve 
size  to  plastic  zone  size.  The  plastic  zone  size  rp  was  estimated  by  measuring  the  extent 
of  the  visual  evidence  of  plastic  deformation  seen  in  the  photographs.  The  initial  curve 
size  r0  was  estimated  using  Eq  (9)  when  r0/rp<\  and  Eq  (5)  when  r0/rp>  1.  It  is  seen  from 
Fig.  14  that  for  r0/rp< 0.35  the  caustics  agree  in  shape  with  the  caustic  predicted  from  the 
plane  stress  HRR  field  (Fig.  2a).  For  r0/rp=  1.4  (Fig.  14f)  the  caustic  has  the  shape  predicted 
using  the  singular  elastic  field. 

Comparing  the  experimentally  observed  caustics  with  the  numerically  generated  ones  of 
Fig.  10  it  is  seen  that  in  both  cases  there  is  a  transition  from  an  "HRR  caustic"  to  an  "elastic 
caustic"  as  r0/rp  goes  from  0.19  to  1.4.  The  transition  away  from  the  HRR  caustic  occurs  at 
ro/rp~0. 35  for  the  experimental  results  and  at  ro/rp~030  for  the  numerical  results.  Thus  it 
appears  that  the  transition  from  the  HRR  caustic  to  the  elastic  caustic  occurs  sooner  in  the 
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numerical  model  than  in  the  experiments.  However,  the  general  trend  is  similar  in  each 
case. 

By  comparing  the  shapes  of  the  experimentally  observed  caustics  to  the  predicted  caus¬ 
tic  shape  (Fig.  2a)  further  information  on  the  region  of  dominance  of  the  HRR  field  can  be 
obtained.  If  the  caustic  is  not  shaped  similar  to  the  predicted  caustic  then  the  displace¬ 
ment  fields  near  the  initial  curve  cannot  correspond  to  that  of  the  HRR  field.  The  caustics 
of  Fig.  14  deviate  from  the  predicted  shape  for  r0/rp> 0.35  indicating  that  under  small 
scale  yielding  an  upper  bound  for  the  dominance  of  the  HRR  field  is  r~0.35rp. 

5.2.  Effect  of  Plasticity  on  Elastic  Caustics. 

Under  conditions  of  plane  stress,  small  scale  yielding  the  elastic  singular  field  dom¬ 
inates  at  some  distance  outside  the  plastic  zone.  Thus  when  the  initial  curve  radius  r0 
satisfies  r0y>rp  and  r0«cj  where  a  is  the  crack  length  or  some  other  relevant  in  plane 
specimen  dimension,  K,  may  be  measured  with  caustics  by  applying  Eq  (4). 

The  results  shown  in  Fig.  8,  show  that  the  elastic  singular  field  does  not  agree  with  the 
actual  stress  field  when  x<1.5rp.  Thus  values  of  K,  measured  with  caustics  are  expected  to 
be  affected  by  the  plastic  zone  when  r0<l.5rp. 

To  quantify  the  effect  of  plasticity  on  caustics  from  outside  the  plastic  zone,  experi¬ 
ments  were  performed  measuring  K,  from  caustics,  Kcaus,  at  different  distances  from  the 
crack  tip  and  comparing  Kcaus  to  K,  measured  from  the  boundary  conditions,  KBC. 

The  results  are  presented  in  Fig.  15  where  Kcails/KBC  is  plotted  vs.  r0/rp.  The  initial 
curve  radius  r0  is  calculated  using  Eq  (5).  For  all  of  the  tests  rp  was  calculated  using  the 
numerical  results 


rp  =  0.25 


Ki 

Oq 


(11) 


for  the  4340  steel  (n=8.7)  and  using  the  estimation 
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(12) 


for  the  1018  steel  (n— ►oo). 

According  to  [15],  for  elastic  solids  plane  stress  conditions  prevail  for  r>0.5h,  where  h 
is  the  specimen  thickness.  Therefore  only  caustics  satisfying  r0>0.5h  were  plotted,  ensur¬ 
ing  plane  stress  conditions.  The  parameter  rp/w  in  the  figure  is  the  ratio  of  plastic  zone 
size  to  uncracked  ligament  size,  w.  This  parameter  indicates  the  extent  of  yielding  that 
has  occured.  No  data  were  obtained  for  rp/w>.046  because  beyond  that  level  of  yielding  it 
was  not  possible  to  observe  a  caustic  that  was  shaped  like  the  predicted  caustic  since 
small  scale  yielding  conditions  were  no  longer  satisfied.  Some  data  points  are  given  for 
0.7<ro/rp<l.0.  It  was  observed  that  for  r0/rp  as  small  as  0.7  the  caustic  curve  remained 
close  to  the  elastic  shape.  Thus  it  is  appropriate  to  present  data  for  ro/rp~0.7  along  with 
data  for  r0>rp. 

The  results  of  Fig.  15  show  that  away  from  the  plastic  zone  Eq  (4)  is  valid.  However,  for 
r0<l.5rp,  Kcaus/KBC  deviates  from  1.0  indicating  that  serious  errors  will  occur  if  Eq  (4)  is 
applied  for  the  evaluation  of  caustics  in  this  range.  This  limit  does  not  appear  to  depend 
on  hardening  for  n> 9  since  the  4340  steel  and  the  1018  steel  data  are  consistent  with  each 
other. 

As  stated  above  the  caustics  retain  the  shape  predicted  by  the  elastic  analysis  even  for 
r0/rp  as  small  as  1.0.  Thus  the  effect  of  the  plastic  zone  on  the  caustic  cannot  be  judged 
by  observation  of  the  caustic  shape.  The  invarience  of  the  shape  of  the  caustic  is 
explained  by  the  numerical  results.  In  Fig.  16  the  angular  distribution  of  the  sum  (Tu+a22  is 
given  for  different  distances  from  the  crack  tip.  It  is  seen  that  an+a22  generally  follows 
the  singular  elastic  field  even  for  r0/rp  as  small  as  1.2.  This  result  is  important  because  it 

is  the  out  of  plane  displacement  u3  = |eru+ <r22j  that  is  relevant  to  caustics,  and 

because  the  shape  of  the  caustics  is  mainly  determined  by  the  angular  distribution  of  u3. 
Since  the  angular  distribution  of  u3  is  in  good  agreement  with  the  singular  elastic  field,  the 
caustic  shapes  do  not  deviate  much  from  the  predicted  shape  as  r0-*rp.  It  is  interesting  to 
note  that  for  small  values  of  r0/rp  the  angular  variation  of  individual  stress  components 
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<7u,a22  does  not  follow  the  elastic  K ,  field  distribution  as  well  as  the  sum  crn+cr22  does. 

5.3.  Three  Dimensional  Effects. 

It  has  been  demonstrated  above  that  the  plane  stress,  HRR  field  dominates  for  r<0.30rp. 
Thus  the  application  of  Eq.  7  for  analyzing  caustics  is  restricted  to  r0/rp< 0.30  if  plane 
stress  conditions  are  satisfied.  However,  due  to  the  finite  thickness  of  an  actual  test  speci¬ 
men  there  is  a  region  near  the  crack  tip  in  which  the  deformation  field  is  three  dimen¬ 
sional.  It  was  shown  in  [15]  that  for  elastic  materials  this  region  extends  for  distances 
from  the  crack  less  than  one  half  of  the  specimen  thickness.  Here  the  three  dimensional 
effects  in  elastic-plastic  materials  are  investigated. 

Experiments  were  performed  measuring  Jcaus  using  caustics  at  different  distances  from 
the  crack  tip  and  then  comparing  Jcaus  to  JBC,  measured  from  the  boundary  conditions.  The 
specimens  used  for  these  tests  are  the  same  specimens  that  were  used  for  investigating 
the  caustics  from  outside  the  plastic  zone.  The  only  difference  is  that  here  the  initial  curve 
was  placed  well  inside  the  plastic  zone,  and  the  plastic  zone  size  was  generally  larger. 

The  collective  results  are  presented  in  Fig.  17  where  Jcaus/JBC  is  plotted  vs.  r0/h.  The  ini¬ 
tial  curve  size,  r0,  was  calculated  using  Eq  (9).  Only  data  points  corresponding  to  caustics 
for  which  the  aspect  ratio  Dy/Dx  was  in  good  agreement  with  the  value  predicted  by  the 
HRR  field  were  plotted.  Here  Dy  and  Dx  are  the  sizes  of  the  caustic  in  the  x2  and  Xj  direc¬ 
tions.  This  was  done  to  insure  that  only  the  three  dimensional  effects  were  being  investi¬ 
gated.  Other  effects,  such  as  the  variation  of  caustic  shape  as  r0/rp  increases,  are  thus 
minimized.  The  maximum  ratio  of  initial  curve  to  plastic  zone  size  was  rQ/rp  <  0.45. 

No  results  are  given  in  Fig.  17  for  r0/h  >  0.6.  The  reason  for  this  is  that  in  a  given  test 
(h,  kp  constant), as  r0/h  is  increased,  r0/rp  is  also  increased.  When  r0/rp  becomes  large  the 
caustic  shape- is  no  longer  similar  to  the  HRR  shape  and  thus  a  calculation  of  J  using  the 
analysis  of  caustics  based  on  the  HRR  field  is  meaningless  in  such  cases.  It  was  found  that 
rp  could  not  be  increased  indefinitely  for  two  reasons.  First,  the  thin  specimens  used  here 
are  prone  to  out  of  plane  bending  and  buckling  at  the  high  loads  needed  to  make  rp  large. 
Second,  extensive  plastic  deformation  reduces  the  reflectivity  of  the  specimens  reducing 
the  definition  of  the  resulting  caustic. 
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The  results  show  that  for  r0/h  <  0.6,  JcaUs/JBc  is  less  than  one,  indicating  that  three 
dimensional  effects  are  important  in  that  region.  The  extent  of  the  region  of  three  dimen¬ 
sionality  seems  to  be  larger  than  the  extent  for  elastic  materials.  The  plane  stress  crack  tip 
solution,  Eq  (3)  predicts  that  u3~* oo  as  r-+0.  However,  since  this  is  impossible,  the  actual 
u3  is  less  than  predicted.  Thus  so  are  the  values  of  JcaUs/JBc  since  the  size  of  the  caustic 
increases  with  increasing  u3  and  du3/dxa. 

The  results  do  not  seem  to  depend  on  specimen  configuration  or  on  the  particular  steel. 
The  reason  for  this  consistency  is  that  the  tests  were  all  carried  out  under  nominally  small 
scale  yielding  conditions,  thus  the  only  in  plane  distance  that  is  important  is  the  plastic 
zone  size.  Also  note  that  the  4340  steel  has  a  relatively  low  hardening  making  it  nearly 
elastic-perfectly  plastic  like  the  1018  steel. 


6.  CONCLUSIONS 

Application  of  the  method  of  caustics  for  the  measurement  of  the  J  integral  in  plasti¬ 
cally  deforming  materials  was  investigated.  It  was  found  that  under  small  scale  yielding 
conditions  the  plane  stress,  HRR  field  dominates  for  r<0.3rp.  Thus  the  application  of  caus¬ 
tics  as  based  on  the  HRR  field  is  restricted  to  initial  curves  satisfying  ro<0.3rp.  For  caustics 
obtained  from  the  elastic  region  the  presence  of  a  plastic  zone  introduces  errors  in  the 
measurement  of  the  elastic  stress  intensity  factor  when  r0<l.5rp.  An  analysis  of  caustics, 
based  on  the  results  of  the  plane  stress,  full  field  numerical  solution,  was  introduced.  This 
analysis  is  valid  for  all  initial  curves.  In  an  actual  test  specimen  there  is  a  region  near  the 
crack  tip  where  the  deformation  field  is  not  one  of  plane  stress  but  is  three  dimensional. 
It  was  found  that  the  extent  of  this  region  for  plastically  deforming  materials  extends  to  at 
least  r<0.6h,  where  h  is  the  specimen  thickness. 
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Figure  1.  Formation  of  caustic  due  to  reflection  of  light  from  the  polished,  deformed 
specimen  surface  near  the  crack  tip. 


Figure  2.  Caustics  formed  due  to  reflection  of  light  from  within  the  plane-stress  HRR 
dominant  zone,  (a)  Predicted  (from  [3])  (b)  Experimental. 
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Figure  3.  Initial  curve  and  plastic  zone  geometries. 
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Figure  4.  (a)  Finite  element  mesh  for  numerical  model,  (b)  Detail  of  mesh  near  crack  tip. 


Figure  5.  Plane  stress  plastic  zone  for  n=9.  (a)  Numerical  (b)  Experimental. 
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Figure  6.  Stress  distribution  ahead  of  the  crack  tip.  r0  =  yield  stress  in  shear. 
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Figure  7.  Plastic  strain  distribution  ahead  of  the  crack  tip. 
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Figure  8.  Stress  distribution  ahead  of  the  crack  tip  compared  to  stresses  from  the  linear 
elastic  Kt  field. 
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Figure  9.  Smoothed  out-of-plane  displacement  field 


Figure  10.  Caustics  generated  by  using  the  out  of  plane  displacements  calculated  from  the 
numerical  model. 
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Figure  1 1.  Relation  between  caustic  diameter  and  J  integral. 


Figure  12.  Experimental  setup  for  photographing  reflected  caustics. 
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Figure  13.  Specimen  dimensions  and  geometries.  AH  dimensions  in  cm. 


Figure  14.  Sequence  of  caustics  obtained  for  reflection  of  light  from  regions  near  a  plasti¬ 
cally  deforming  crack  tip. 
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Figure  15.  Kcaus/KBC  vs.  r0/rp.  Deviation  from  1  0  indicates  error  caused  by  plasticity. 


Figure  16.  Angular  distribution  of  (un+(T22)  for  different  distances  from  the  crack  tip. 
Solid  line  is  stress  distribution  of  K[  field. 
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Figure  17.  Jeans  I  he  vs.  r0/fi.  Deviation  from  1.0  indicates  error  due  to  three  dimensional¬ 
ity  for  r0/h< 0.6. 


